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Abstract 

In the context of gauge/gravity duals with flavor, we examine heavy-light 
mesons which involve a heavy and a light quark. For this purpose we embed two 
D7 brane probes at different positions into the gravity background. We establish 
the non-Abelian Dirac-Born-Infeld (DBI) action for these probes, in which the 
U(2) matrix describing the embedding is diagonal. The fluctuations of the brane 
probes correspond to the mesons. In particular, the off-diagonal elements of the 
{7(2) fluctuation matrix correspond to the heavy-light mesons, while the diagonal 
elements correspond to the light-light and heavy-heavy mesons, respectively. The 
heavy-light mesons scale differently with the 't Hooft coupling than the mesons 
involving quarks of equal mass. The model describes both scalar and vector 
mesons. For different dilaton-deformed gravity backgrounds, we also calculate 
the Wilson loop energy, and compare with the meson masses. 
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1 Introduction 



Recently, based on the gauge/gravity correspondence [T], many non-perturbative prop- 
erties of Yang-Mills theories with quarks have been uncovered in terms of superstring 
theory [2]-[TT]. Flavor quarks in the fundamental representation of the gauge group are 
introduced by embedding one or several probe branes into an appropriate bulk gravity 
background, in order to describe large N gauge theories similar to QCD. Many suc- 
cessful results have been obtained for the properties of quarks and their bound states: 
Mesons spectra have been studied by many authors, including [3]-|20j. The mass spec- 
tra of fermionic operators with fundamental fields ( "mesinos" ) have been discussed in 

[223- 

Up to now, most of these investigations have been devoted to the case that the D7 
branes are embedded at the same place, such that the flavor group forms U(Nf) for Nf 
D7 branes. Effects of the non-abelian nature of U(Nf) for Nf > 1 in relation to bulk 
instantons have been studied in [23J and in [24j-[26j. Examples of D7 branes embedded 
at different positions occur in studies of meson decay via string breaking [27]-[29]. 

Here we propose a holographic model for heavy-light mesons based on a non-abelian 
Dirac-Born-Infeld (DBI) action. In this model, two D7 brane probes are embedded at 
different positions, such that they provide different quark mass states depending on 
the flavor. It is possible to choose any mass difference by separating the branes by an 
appropriate distance. In terms of this model, we study meson states with heavy and 
light quarks, such as mesons with charm or bottom degrees of freedom for instance. 
Heavy-light mesons with large spin have been studied in (3TH I3T] . 

In [32] , an effective model for heavy-light mesons such as the B meson has been pro- 
posed which is based on the Polyakov action. This model gives a qualitative description 
of the B and excited B* states. It also gives rise to a dependence of the heavy-light me- 
son mass on the 't Hooft coupling (A) of the form Mhl/^h = 1 + const/vA + C(A _1 ), 
with Mhl the heavy-light meson mass and m# the heavy quark mass (The light quark 
mass m L has been set to zero here.). This differs from the 't Hooft coupling dependence 
of the light-light or heavy-heavy mesons [3J, for which M ex. m/ y/\. However the result 
of [32] ensures that the heavy-light meson mass equals the heavy quark mass in the 
large A limit, Mhl — n^H, as expected from observation and dimensional analysis. 

Here we present an alternative approach to the same heavy-light meson mass in 
terms of a non-Abelian Dirac-Born-Infeld (DBI) action for two D7 branes at different 
positions. We use the non-Abelian DBI action for curved lOd backgrounds which has 
been proposed by Myers [33]. In this action, the world- volume fields are assigned to 
U(Nf) matrix-valued functions for Nf D7 branes. We choose Nf = 2. The embedding 
configuration of the two D7 branes is determined by the diagonal components of the 
scalar fields. The corresponding equation of motion is solved by the profile functions of 
two separated branes, one of which corresponds to the heavy and one to the light quark. 
The quark masses are given by the boundary values of the two embedded branes. The 
fluctuations of the diagonal elements of the 2x2 flavor matrices correspond to the 
light-light and heavy-heavy mesons, respectively. On the other hand, the off-diagonal 
components of the fluctuations of the fields on the branes are identified with the heavy- 
light mesons. We calculate the spectrum of these heavy-light mesons. 

For the A dependence of the heavy-light meson mass, we find that it is similar to the 
one observed using the Polyakov action approach in [32]. A finite contribution to the 
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mass remains in the limit of A — > oo. This contribution corresponds to the minimum 
energy of a classical string connecting two separated D7 branes, and thus is equivalent 
to the mass obtained from the Polyakov action. 

This A dependence persists if we consider the D3+D(— 1) gravity background of [33] . 
In the field theory dual to this background, a condensate q = n 2 (F 2 ) is switched on. 
D7 embeddings and chiral symmetry breaking for a non-supersymmetric version of this 
background have been studied in [10]. The A dependence is very similar to that in the 
supersymmetric background. 

It is instructive to compare the A dependence of the meson spectra with the A 
dependence of the tension. For a classical string stretched between the two D7 brane 
probes, the string tension is independent of A, in agreement with the heavy-light meson 
mass result found in [32] as well as in the present paper. For heavy-light mesons this 
tension contributes to the meson mass even if the distance L between the quark and 
anti-quark in the four-dimensional boundary space is zero, in which case it contributes 
E = rax — rriL to the Wilson line energy. For the heavy- heavy and light-light mesons, 
the string tension scales as m; ]/y/X for small L [31 HD]- At large L, when the dual gauge 
theory is in the quark confinement phase, there is a long range linear potential for 
all the mesons considered. For the heavy-heavy and light-light mesons this was found 
in |10]. 

Our model allows to describe both scalar and vector mesons. In this respect, it goes 
beyond the effective model of [32]- For the supersymmetric D3 + D(— 1) background 
of [31], we find that the vector and scalar meson masses differ, due to the different 
dependence of the fluctuation equations on the non-trivial dilaton. Note that when 
adding D7 probes to the D3 + D(— 1) background, supersymmetry is broken to Af = 1 
and the vector and scalar mesons are not in the same multiplet any more. For very 
large heavy quark mass, M = 2 supersymmetry is restored, and the vector and scalar 
meson masses become degenerate again. This is consistent with the phenomenological 
fact from heavy-quark theory that spin effects are suppressed by powers of the inverse 
heavy quark mass. Of course, here this is due to M = 2 supersymmetry restoration. 
We leave an investigation of this mechanism for non-supersymmetric backgrounds to 
the future. 

Our paper is organized as follows. In section 2, we give the non-Abelian DBI action 
and the D7 brane embedding model is proposed. In section 3, the number of D7 
branes is restricted to two (Nf = 2) and the action is expanded by fluctuations to see 
the meson spectra, which are shown in the section 4 for the case of HL mesons. In 
section 5, the potential between quark and anti-quark is given through Wilson loop. 
The summary is given in the final section. 
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2 Embedding of D7 branes 



2.1 Non-Abelian Dirac-Born-Infeld action 

We start from the non-Abelian Dirac-Born-Infeld action proposed by Myers in |33j . 
This action describes the dynamics of Nf Dp-branes in a background with metric G mn 
and is given by 

S Nf = -r p J tP+^e^STr (V- det(P[G rs + G ra {Q-^ - 5)" b G sb ] + T^F rs ) ^/drtQ^j , 

(1) 

where the matrix Q a b is defined by 

Q a b = 5 a b + tT[X a ,X c }G cb (2) 

where T _1 = 27ra', and X a are the coordinates transverse to the stack of branes, which 
now take values in a U(Nf) algebra. The symbol STr denotes the symmetrized trace 
STr(Ai...A n ) = ^Tr(Ai...A n + all permutations) and is needed to avoid the ambiguity 
of the ordering of the expansion of all fields in the DBI action. [35J. 

In our convention, r, s — 0, 1, ...,p and a, b — p + 1, 9 label the world- volume di- 
rections and the directions transverse to the Dp-branes, respectively; m, n = 0, 1, • • • ,9 
are the lOd spacetime indices. P[a rs ] denotes the pull-back of a lOd tensor a mn to the 
world-volume of the branes. A peculiarity of the non-Abelian DBI action is that the 
pull-back matrix is given by the covariant derivative 

D r X a = d r X a + i[A r ,X a ] , (3) 

with partial derivatives d r = d/d£, r , non-Abelian world- volume gauge field A r and 
transverse coordinates X a . F rs is the corresponding world- volume field strength. 

For diagonal brane embeddings the commutator [X a ,X 6 ] is small, as can be seen 
as follows. Set X a as X a = X a + a , where a denotes small quantum fluctuations 
around a diagonal embedding matrix X a . Then we have 

[X\ X b ] = [X a , X b ] + [X a , <p b ] + [<p\ X b ] + [<p a , <p b ] , (4) 

where the first term [A a ,X' ) ] vanishes and the remaining terms are small. Next, as- 
suming also a diagonal metric G mn and employing the approximation (Q -1 — S) ab ~ 
— iT[X a , X b ], we rewrite the pull-back in the action ([1]) as 

P[G rs + G^Q- 1 - 5) ab G sb ] « G rs + D r X a D s X b (G ab - iT[X c , X d ]G ac G bd ) . (5) 

Then, the action (jTJ is expanded in powers of [X a ,X b ] up to 0(X 4 ) as 

S Nf = r p J al p+1 i e^STr | y/- det(G rs + G ab D r X a D s X b + T- l F rs ) 

x (l-\(TG ac [X^X b }) 2 ^ } . (6) 

The factor in the second line descends from the expansion of \/det Q a b . For a flat 
spacetime background G mn = r] mn , the action ([6]) agrees with that found in [35J. 
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2.2 Diagonal embedding ansatz in non-Abelian DBI action 



The non-Abelian DBI action will now be used to find the embedding of Nf probe D7 
branes in different gravity backgrounds. The embedding profiles correspond to the 
classical solutions for the scalar fields in the D7 brane action. In our case, the scalar 
fields X a are U(Nf) matrix valued functions which makes it difficult to obtain a general 
form of the profile functions. In order to simplify the problem, we use the diagonal 
ansatz 

X = diag«--- ,w a Nf ), (7) 

thereby setting all off-diagonal components to zero. Here each of the functions wf 
corresponds to one of the Nf D7 branes. - It would be an interesting problem to also 
include the off-diagonal components of X a and to solve the corresponding embedding 
equations. For example, for a non-trivial world-volume gauge field F rs the off-diagonal 
components provide a Bl-on configuration which connects two branes [3S]- We postpone 
the discussion of such configurations to the future. 

The quark mass for each flavor is given by the asymptotic value of wf in the ultra- 
violet limit. They are the integration constants and given by hand as parameters of 
the theoryEl The equations of motion for the wf are obtained from the action 

S Nf = r 7 j d 8 £ e~*STr (^/ - det(G rs + G ab d r wfd s w\)^ 

N f 

= r 7 f d 8 £ e"* \J-^{G rs + G ab d r wfd s wt) (8) 
i=i 

which is Eq. (Ej) for the embedding (JTj) and p = 7. The essential point is here that for 
the diagonal ansatz (JTj), we obtain Nf decoupled equations of motion for the wf such 
that the embeddings of each of the probe branes is independent of the other. In other 
words, for diagonal embeddings the non-Abelian DBI action reduces to the sum of Nf 
abelian DBI actions[j] 

If there are additional fluxes in the background, the action (jSj) must be supple- 
mented by appropriate Wess-Zumino terms. 



2.3 Specific gravity background 

Let us now find the embedding functions w = Wi (i = 1, Nf) of the probe branes for 
some phenomenologically interesting supergravity backgrounds. 

As a specific supergravity background, we consider the following lOd background 
in string frame given by a non-trivial dilaton $ and axion \ [34, 10J, 

ds 2 10 = e* /2 (f-A^r^dxW + ^dr 2 + R 2 dn\\ . (9) 



*The asymptotic values of (possible) off-diagonal components, say wfj, represent the mass-mixing 
of different quark flavors i and j. We neglect these here. 

^This is only true, if we ignore the fluctuations around the brane embeddings. As we will see in 
Sec. 3, the non-Abelian DBI action also describes fluctuations of strings stretched in between two 
different branes. 
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Here two typical solutions are considered. One is the supersymmetric solution 



A — l, e* = 1 + 



Q 



X 



-e + Xo 



(10) 



and the other is non-supersymmetric and given by 



A(r) = (l - ft 8 ) , X = 



mi 



>Ao) 4 -i, 

The first solution is dual to M = 2 super Yang-Mills theory with gauge condensate q 
and is chirally symmetric. Both supersymmetry and chiral symmetry are broken for 
the second solution [TU1. 



In order to obtain the induced metric on the D7 world-volume, we rewrite the 
six-dimensional part of the metric ([9]) in the form 



^dr 2 + R 2 dn 2 5 = % {dp 2 + P 2 dn 2 3 + (dx 8 ) 2 + (dx 9 ) 2 ) 



(12) 



where r 2 = p 2 + (X 8 ) 2 + (X 9 ) 2 . Due to the rotational invariance in the X s — X 9 plane, 
we may set X 8 = and X 9 = w(p) without loss of generality. Then the induced metric 
on the D7 brane is given by 



ds\ = e*/ 2 ^A\ u dx^+ ^ ((1 + (d p w) 2 )dp 2 + p 2 dnj) | . 



(13) 



In static gauge the action for the D7 probe is given by 

s D7 = s DBI + s wz = -r 7 j d 8 z^r 3P 3 (VeVi + K) 2 - C *) ' ( 14 ) 

where C$ = q/r A = — 1 denotes the Wess-Zumino term coming from A$, the Hodge 
dual of the axion [TO] . The Wess-Zumino term is only required for the supersymmetric 
background ffTUl) . The equation of motion for the embedding function w(p) is 



w 



p + w w 

1 

4 



7V / TTKj 2 ($ + 4iogA) / 



for the non-supersymmetric case (fTTj) . and 



w'[ - + ($ + 41ogA)'J + i 



w 



/^2 



w 



p + w w' 



1 - v 7 ! + K) 2 



l + (w 



l\2 



p / 1 + (w') 2 



(15) 



(16) 



for the supersymmetric case (jlOl) . Here the prime denotes the derivative with respect 
to p. 

In deriving these equations, we have to take into account that r = \/ p 2 + w(p) 2 . 
We therefore have to extract the variation of w also from the functions of A(r) and 
$(r). For example, the variation of A(r) with respect to w is obtained as 



SA(r) 



Or 2 
dw 



d r 2A(r)5w + 



w 



p + wd p w 



d A5w + ■ ■ 



Here, the expression after the second equality sign shows the change of variable from r 
to p in the derivative. The prefactor of the first term of (fl5|) and (fl6l) originates from 
this variable changing procedure. 

Solving the above equation for w, we find the profile functions of the D7 brane 
embedded at the separated places and then we find simultaneously the quark properties, 
the quark mass m q and the chiral condensate (\E , \I / ), where \l/ denotes the quark field. 
The details of the solutions are shown in [TO] . 



3 D7-brane fluctuations 

In this section we derive the actions of the scalar and vector D7 brane fluctuations dual 
to heavy-heavy, light-light and heavy-light mesons. These actions will be used in the 
next section to find the fluctuation spectrum in the backgrounds of the type (J9j). 

For this, we return to the brane action ([6]). At this stage, we restrict to the case 
of Nf = 2 flavors or two D7 branes such that the scalar and vector fields in the non- 
Abelian DBI action are represented by 2 x 2-matrices. For the classical embedding, we 
choose the diagonal configuration given by 

X 8 = 0, X 9 = ( 7 ° V (17) 
\ w 2 J 

In terms of the Pauli matrices 



2 V 1 7 ' ' 2 V 1 ) ' ' 2 V i J ' ' 2 V -1 

X 9 can be rewritten as 

X 9 = wt + vt 3 , w 1 = (w + v)/2, w 2 = (w-v)/2, (19) 

where v = Wi — w 2 . The asymptotic boundary values of w\ and w 2 correspond to the 
heavy and light quark masses, respectively. When v = 0, the two branes are at the 
same place, W\ = w 2 = w, corresponding to a £7(2) flavor symmetry. For v ^ this 
flavor symmetry is explicitly broken. 

The scalar and gauge field fluctuations are taken to be of the form (a = 8, 9) 

X 9 = X 9 + (j) 9 , X 8 = 8 , (20) 

a = 0^° + 2 V, A r = A r T° + Ay, (21) 



and can be written as 



%i <f>° 



(22) 



similarly A r . The diagonal elements <j)± = 0q ± 0| describe fluctuations of each brane 
and are dual to the heavy-heavy and light-light mesons. On the other hand, the off- 
diagonal elements 4>\ 2 = (p\ — icfi?, and 4>2i = 4>1 + correspond to fluctuations of 
strings stretched between the two branes and are dual to the heavy-light mesons. The 
mass of this last type of fluctuations will depend on v. - A similar structure emerges 
also for gauge field fluctuations A r , as discussed below. 
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These meson mass spectra are obtained by solving the linearized equation of motions 
for the field fluctuations. Using the expansions (120]) and assuming small fluctuations 
<J) a and A r , the action (ED is rewritten as 



S Nf =2 = r 7 J d s £ STr |e" V~ det(a rs ) + G 88 G 99 ^ ((0 8 ) 2 + (0 8 ) 2 ) - 2 



^ , (23) 



where 

a rs = G rs + G a bD r X a D s X b + T 1 F rs . (24) 

and 4>\ and 0| as defined in (12T1) . The other components of 8 and 9 do not appear 
explicitly, but contribute to a rs in (123]) . 

When evaluating (|23|) . we have to take into account that the radial coordinate r, 
which occurs in G rs , G a b and $ on the brane, is a matrix- valued function. Including 
the fluctuations, r is of the form 

r 2 = p 2 + (X 8 ) 2 + (X 9 ) 2 = p 2 + ( W r + t;r 3 + 9 ) 2 + (0 8 ) 2 (25) 

and, for instance, G rs denotes the matrix 

G rs ^ ( ( l rs \ r=ri1 Jj-j™ ) , (26) 

where (i, j = 1, 2) are the matrix elements of r. 
It is convenient to write a rs in (1231) as 



where 

d rs = G r s G 99 r X 9 s X 9 . (28) 

The explicit fluctuation dependent part 5a rs is expanded in terms of the power series 
of the fluctuations as 

Sa rs = a® +4 2 J + (29) 

where 

a« = G 99 (d r X 9 d s <f) 9 + d r <f) 9 d s X 9 - i[A r , X 9 ]d s X 9 - id r X 9 [A s , X 9 ]) + T^F^ , (30) 
ag> = G 99 (-z[A r ,X 9 ]9 s 9 - td r <J) 9 [A s ,X 9 } + d^d^ 9 - [A r ,X 9 ][A s ,X 9 ] 

-i[Ar, <P 9 }d s X 9 - id r X 9 [A s , 9 ]) + G 88 <9 r 8 «9 s 8 . (31) 

Moreover, since the metric depends on r given by (125!) . a rs in (|28|) still includes (im- 
plicitly) the matrix- valued quantum fluctuations 8 and 9 . 

Then the action §6§ is expanded up to quadratic order in the fluctuations, 

S Nf=2 = r 7 J d 8 £ STr je-V-det(a rfl ) + ^(a" 1 ^) + ^ (tr^a" 1 ^)) 2 
-itr rs ((rt«) 2 ) + ^tr rs (a-V 2 )) -^G 88 G 99 ((0 8 ) 2 + (^)> 2 + ■ ■ ■ ) ) . (32) 
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Lagrangian for scalar fluctuations 

From (13TI) . we see that the fluctuations <fi 9 and A r are mixing, while 8 does not mix with 
any other field. For simplicity, we consider only the 8 fluctuations. When evaluating 
the symmetrized trace in (I3"2j) it has to be kept in mind that the diagonal flavor matrix 
elements of the embedding have to be evaluated at W\(p) and w 2 {p), respectively. To 
quadratic order the Lagrangian for the 8 fluctuations reads 

4? = \d r ,F\ m ((0 8 ) 2 + (0 8 ) 2 + (0 8 ) 2 ) + \d r >F\ W2 ((0 8 _ f + (0 8 ) 2 + (0 8 ) 2 ) 
+ ~(F G 88 a rs )\ Wl d r <f> 8 + d s <f>l + ~(F G m a rs )\ W2 d r (jtd s 4t 

o o 

+ l - ((F G 88 a rs )\ Wl + (F G 88 a rs )\ W2 ) (<^<9 S 8 + d r <f 2 d s( f>l) 

- V - ((F G 88 G 99 )\ Wl + (F G 88 G 99 )\ W2 ) ((0 8 ) 2 + (0 8 ) 2 ) , (33) 

o 

where 

F = e^y / -det a rs . (34) 

In the above equation (|33j) . the notation K(r)\ w . means that r in any function K(r) 
is replaced by f{ = \J p 2 + wf, K(r)\ Wi = K(^J p 2 + wf). Moreover, the first two terms 
of 0331) are derived from the matrix- valued coordinate r, which includes the fluctuations 
(f) s and (f) 9 . The first two terms of (1331) are essential for finding the correct spectrum. 
In particular, in the non-supersymmetric case these terms are crucial for finding the 
Nambu-Goldstone bosons of chiral symmetry breaking. 

From the Lagrangian (133]) we obtain the standard equations of motion for the 
heavy-heavy and light-light modes (f)\ as well as a new equation for the heavy-light 
fluctuations. Note that the mass spectra of 8 2 depend on both profile functions 
^i,2(p)- ~ For the details of the meson spectrum in the single brane case see [37] . 

Lagrangian for vector fluctuations 

For the vector meson, we find a similar separation of the modes as for the scalars. 
However there is a mixing of the A].' 2 and 9 modes, as we now discuss. For the case of 
constant w 1 , we can see that this mixing can be removed by the gauge transformation 
of Al' 2 . By taking 

Al = Al + -d r <f>l, 
v 

A 2 = A 2 r - l-drfi , (35) 
the mixing part in (!3T!) is written as 

41° = G " (-i[A r ,X 9 ]d s <p 9 - id r 4 9 [A s ,X 9 ] + d^ds^ 9 - 2[A r , X 9 ] [A s , X 9 ]) . 
Again, when explicitly writing out the flavor matrix, we have a diagonal form 

= <*{( m e M ) - ( J S ) T + ^) } ' « 36 > 
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where the top left entry is evaluated at W\ and the bottom right at u> 2 - From f[3"oT) we 
see that the kinetic terms of <fi\ 2 are eliminated by changing A]: 2 to the new variables 
A],' 2 . On the other hand, due to gauge invariance, the kinetic term of A r does not 
give rise to any new kinetic terms for 2 . Instead, new mass terms are generated for 
A). )2 as shown above. The two scalar components 4>\ 2 are gauged away to produce the 
longitudinal component of the vector A^ 2 . 

This is the well-known Higgs mechanism, with X 9 the Higgs scalar and the A®, the 
SU (2) gauge fields. X 8 is not involved in this Higgs mechanism of the gauge symmetry 
breaking. However X 8 is associated with the Nambu-Goldstone mode of the geomet- 
rical U(1)a chiral symmetry breaking. These two symmetry breaking mechanisms are 
not related to each other. 



For constant 101,2, the vector meson part of f[3"2l is given by 




where in flavor space 

af = F rs (A)(l J ) , (38) 

a f = G " ( I °i ) T (% A > + A * A *) ■ (39) 

Then we obtain 

+ 1 ((FcW)U + (Fa rp a^)\ W2 ) + F 2 S F^) 

2 

+ j ((F G 99 a rs )\ wi + (F G 99 a rs )\ W2 ) (ijij + A 2 r A 2 )j , (40) 

where 

F ± = F° ± F 3 . (41) 

rs rs rs \ J 

Again the dependence on the position of the two distinct branes, characterized by the 
two profile functions Wi(p), i — 1, 2, is stated explicitly as in fl33l above. 

We give a brief summary of this section, (i) For heavy-heavy or light-light mesons, 
the modes of 0q 3 and A r QZ recombine into (f>± and A±, and their spectra coincide with 
the one obtained from one individual probe brane. (ii) For heavy-light mesons, only 8 2 
and A\ 2 remain as independent fluctuation variables. (f)\ 2 are gauged away and A\ 2 
become massive as a result of the Higgs mechanism, (iii) The mesons corresponding to 
(f)\ 2 have a mass which depends on both brane embeddings w± and u> 2 . These modes 
are not affected by the above Higgs mechanism. 
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4 Mass spectra of heavy-light mesons 



In this section we examine the mass spectra of heavy-light (HL) mesons in the two 
background solutions. 

4.1 Supersymmetric case 

First, we consider the supersymmetric background solution given by ffTUj) . 
4.1.1 Scalar fluctuations 

For the supersymmetric solution ffTUj) . the equation of motion for the fluctuations (f)\ 2 
dual to HL mesons is obtained from (|33|) . We find 



ag + Vjg+2 + + (42) 

p 9 p 2 e*( ri ) + e*( r2 ) ' 



Mi = e *<"> f 5 p $(rOS P + ^ i? 4 J , (43) 

where rf = p 2 +w 2 (i = 1, 2) and 2 are denoted by 0. Here — /(/ + 2) is the eigenvalue 
of the Laplace operator on S 3 . 

In order to study the qualitative behavior of the spectrum, we first consider the 
non-confining case q = 0. In this case, the dilaton is trivial, $ = 0, and Eq. 
simplifies to 




Q2 3 a l{l + 2) M 2 -v 2 ( R 2 \ 2 ( R \ i i , 

^V'-V^^IAtt^J Htt^i' D 0=o - (44) 

The mass of the fluctuations M is interpreted as the heavy-light meson mass in the 
dual gauge theory. Since w\ and w 2 are mixed in a complicated way, we must solve 
this equation numerically. 

Before proceeding to the numerical solution, we consider two special cases in which 
the spectrum can be determined analytically. For w\ = w 2 = w, we get v = and 
the equation reduces to the one given by Kruczenski et al [i3J which can be solved 
analytically. In this case, the meson masses can be expressed in terms of the quark 
mass m = w/(2Tca f ) and the 't Hooft coupling A = R 4 /47ia' 2 by 

2 

777 

M 2 = 4tt— — (n + I + l)(n + I + 2) , (45) 
A 

where n denotes the node number of the eigen functions for I = and and I represents 
the angular momentum of 5" 3 in the world-volume of the D7 brane. This spectrum 
represents therefore the one for the heavy-heavy and light-light mesons. 

The other regime in which fj44l) can be solved analytically corresponds to a heavy- 
light meson with a very heavy quark, w 2 W\. In this case the term in (|44p involving 
w 2 is much smaller than the one involving W\ and may be neglected. Eq. (jUj) ap- 
proaches then the equation for a meson with a single flavor [3| in which we replace 

M 2 - v 2 

M 2 — > M 2 = , w w t . (46) 
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1.25 1.5 1.75 2 2.25 2.5 1.5 2 2.5 3 3.5 4 



Fig. 1: Numerical plots of the heavy-light meson mass for different values of the 't Hooft 
coupling A. Here we set T = 1 and = W\ — 1. The dashed curve shows the 
analytical solution (|47l) for A = 3 4 . 



Substituting this into (j4"5|) for n = I = 0, we find 

where we reintroduced the string tension T = l/(27ra') (which was set to one above) 
and defined the quark masses = w 1)2 /(27rQ! / ) as the distances in units of T. 

Eq. ( 1471) implies that the mass of HL mesons has two different contributions. The 
first term proportional to has the same dependence on the 't Hooft coupling as 
in the single flavor case [3]. The second term is dominant at large 't Hooft coupling 
(A — > oo), where the mass of the HL mesons is approximated by the second term, 



2na' 



In this strong-coupling regime, the heavy-light meson mass depends solely on the dif- 
ference of the two quark masses. This is consistent with the result obtained in |32j . 
and provides a lower bound for the HL meson mass. 

Let us now solve Eq. (jHj) numerically and compare the results with the analytical 
solutions. Fig. [1] shows numerical plots of the heavy-light meson mass M^ L (in units 
of A) versus the heavy quark mass m H for various values of the 't Hooft coupling. 
As an example, consider the asymptotics of the graph for A = 3 4 (Fig. [TJ-left). For 
small quark mass differences, m# ~ (ro/, = 1 in the plot), the graph behaves 
very similarly to the heavy- heavy graph (j45p with m = ma- For large quark mass 
differences, m# 3> mi, the numerical graph asymptotes to the analytical curve given 
by Eq. fj47j) (dashed curve). We observe that the approximation ( |4"T|) works particularly 
well for large masses of the heavy quark. 

In Fig. [T]-right we compare the heavy-light meson mass with the heavy-heavy meson 
mass. The essential difference comes again from the 't Hooft coupling dependence. For 
small 't Hooft coupling the HL curve lies below the HH curve. However, there exists 
a critical value of the 't Hooft coupling for which the HL mesons are heavier than 
the HH mesons. In other words, at strong coupling (at large A) and for fixed quark 
masses ma l, the HL meson mass is much larger than the corresponding HH meson 
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mass. This is of course unphysical from the point of view of QCD. However, this seems 
to be a general feature of gravity dual heavy-light models, since the 't Hooft coupling 
dependence found here coincides exactly with the one found in [32J. 

Next, as in [53] . we turn on a gauge condensate q ~ (F^) such that the dual 
supersymmetric gauge theory becomes confining, but remains chirally symmetric. In 
general, the gauge condensate depends on the 't Hooft coupling and we choose q = qXa' 4 
as in [33] . For the dilaton given by ffTUj) , Fig. [2] shows the HL and HH mesons as a 
function of the heavy quark mass for fixed value of the light quark mass mi = 1. We 
observe that the presence of q increases the HL meson masses. This is due to increase 
of the dilaton in the term v 2 e* in the presence of q. The dilaton as given by (TlO|) is 
responsible for quark confinement and an increase in binding energy. Note that the 
term v 2 e* is independent of A. On the other hand, for the HH mesons the q dependence 
via e* disappears at large quark mass m#. This is seen from their equation of motion 
which is obtained from Eq. (1421) by replacing the last term by 



M 1 + M 2 



where 



M 2 

d p $(r)d p + —R 



Aqp 



{w 2 + p 2 ){q+{w 2 + p 2 ) 2 )' ^ 

This shows the g-dependence explicitly. We see that this dependence is small for 
large w, i.e. for large quark mass, for all values of p. Thus in the large w limit, the 
HH meson masses take their M = 2 supersymmetric value of the q = case. This 
implies an increase in the difference between the HL spectrum and the HH spectrum 
for sufficiently large A. This is due to the term v 2 e* which contributes only to the HL 
mass. 



AM 





3.5 



m H 



Fig. 2: Meson masses for non-zero q. The red and blue curves show Mhl, Mhh for 
A = 3 4 , q = 1 (left) and q = 5 (right). The grey curves show the corresponding meson 
masses for q = 0. The presence of q increases the HL meson masses. The lambda 
dependence remains unchanged. 



4.1.2 Vector fluctuations 

From (1401) . we obtain the equation of motion for the vector fields for the confining 
supersymmetric background of [33] • Here we consider the four dimensional vector A 1 ^ 2 , 
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for the case that the components A„ = A4 = are zero. By imposing the gauge 



condition d^A^ 



0, the vector fluctuation equation of motion is 



aft + -a, 



1(1 + 2) Mf + M£ 



+ 



M 2 



2 



A, 



-R A 



r 



(50) 



(51) 



where i — 1,2 and A 1 ^ 2 are denoted by A^. Note the different dilaton dependence as 
compared to the scalar equation ( )44l) . The different dilaton dependence of vector and 
scalar mesons is consistent with the fact that when adding D7 probes to the D3+D(— 1) 
background, supersymmetry is broken to Af = 1. 




1.5 



2.5 



3.5 



Fig. 3: Vector and scalar masses as function of the heavy quark mass for R = 2, q = 10. 
The vector masses are larger than the scalar masses. For large heavy quark mass, they 
become degenerate again. Red: scalar heavy-heavy meson. Red dashed: vector heavy- 
heavy meson. Blue: scalar heavy-light meson. Blue dashed: vector heavy-light meson. 



We find that the vector masses are slightly larger than the scalar masses. For large 
heavy quark masses, the vector spectrum degenerates with the scalar spectrum again, 
since the dilaton dependence becomes negligible. For vanishing dilaton $ = 0, ( 1501) 
reduces to the same form as the scalar (144p . The degeneracy of vector and scalar masses 
in this limit is expected since for very large heavy quark mass, M = 2 supersymmetry 
is restored. This is consistent with the phenomenological fact from heavy-quark theory 
that spin effects are suppressed by powers of the inverse heavy quark mass. Of course, 
here this is due to M = 2 supersymmetry restoration. It would be interesting to 
compare the scalar and vector sectors for a non-supersymmetric gravity background. 
However, since the calculations are much more involved, we leave this for future work. 
In the next section we consider the non-supersymmetric case for just the scalar sector. 

4.2 Non-supersymmetric case and chiral symmetry breaking 

The above analysis is performed for a supersymmetric background, so w is constant 
and there is no chiral condensate. On the other hand, in a background dual to a non- 
supersymmetric theory with chiral symmetry breaking, the profile functions are 
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not constants but vary with p. As a result, the mass spectrum is modified due to the 
presence of the chiral condensate and the related background configuration. 

We now consider the non-supersymmetric background (llip . The corresponding 
mass spectrum of single-flavor mesons has previously been studied in [37J. In the 
following we compute the HL spectrum dual to the scalar fluctuations (p\ 2 and compare 
it with the corresponding HH spectrum. We show that at strong 't Hooft coupling the 
heavy-light meson masses lie below the corresponding heavy-heavy meson masses in 
agreement with phenomenological expectations. 

The linearized equation of motion of the <p — <j>\ 2 fluctuations is given by 



dl + -d p + 



N {1} + N {2} 
F {1} +F {2} 







with 



i% = {F (d p (\ogF)d p + K) — v 2 G} 



G = e™A\ll + {d p wy ^ 



(52) 



(53) 



K = (1 + (d p w) 2 ) - -2K r ) , K r = d r2 log(e*A 4 ) , ( 5 I ) 



(55) 



VI + (^) 2 ' 

and the dilaton $(r) and warp factor A(r) as in (fill . The index {i} (i = 1,2) means 
that we have to substitute either of the profiles Wi^{p)- In the following we consider 
fluctuations with quantum numbers n — I — for the sake of simplicity. 




Fig. 4: Embedding solutions w(p) in the nonsupersymmetric background (TTTT) . 



For the above non-susy configuration, the parameter r plays a role similar to the 
infra-red cut-off Aqcd in QCD. Recall that the background has a singularity at r = r$ 
and is well-defined only for r > r corresponding to energies above Aqcd- Quite 
generally, we expect the parameter r to depend on the (asymptotic) AdS radius R. 
For the computation of the meson spectrum, we will make the simple choice 

r = R. (56) 
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The dependence of the parameter r on R is motivated by the fact that r is also related 
to the gauge field condensate (F 2 ), as can be seen by expanding the dilaton as 

e< t> = l + g/r 4 + • • • , (57) 

where q = v^6^o- Now, q is related to the gauge field condensate (F 2 V ) by q = X(F 2 ), 
see for instance [34] , and thus r oc R. We should notice here that A is running in the 
present case, since gy M = e®^ depends on r. 

Let us now find the numerical fluctuation spectrum. As in the single-flavor case, 
one first computes the D7 brane profiles Wi^p) by numerically solving the embedding 
equation (fl5l) . The asymptotic values of wi^ip) are fixed by the quark masses rriL,H- 
The values of the quark condensate c = (^0) are obtained by requiring the regularity of 
the solutions for all values of p. For the background (jTTj) . the embedding solutions w(p) 
have been found in [TO] and are shown in Fig. H]for various quark masses m [TO]. As 
in [5], the solutions are repelled from the singularity at r = R = 1. Again, fluctuations 
of strings stretching in between two different branes (such as the blue HL string) are 
dual to HL mesons. Strings starting and ending on the same brane (e.g. the red HH 
or the green LL string) correspond to HH or LL mesons. 

The embedding solutions must then be substituted into the equation of motion (152]) . 
Solving this equation for <p, we obtain a numerical spectrum M(m#, m^, R) of HL 
mesons. Fig. shows the resulting HL spectrum in dependence of the heavy quark 
mass run- The light quark mass is set to zero, ttil = 0, and R is kept fixed here (R = 1). 
For run = 0, we recover the (massless) Nambu-Goldstone boson in the spectrum which 
is expected from spontaneous breaking of the £7(1) chiral symmetry [5], [10]. Both the 
HL and the HH spectrum satisfy the Gell-Mann-Oakes-Renner relation M 2 oc m# for 
small m H . 



M 2 




m H for zero light quark mass and R — 1, 2na' = 1. 

We also find that the HL meson masses lie below the HH spectrum, at least for an 
intermediate value of the 't Hooft coupling of 1 < A, as can be seen from Fig. [61 As 
in the supersymmetric case, the heavy-light mesons scale differently with the 't Hooft 
coupling than the heavy-heavy and light-light ones. There exists a critical value for 
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Fig. 6: 't Hooft coupling A dependence (for lua' = 1, i.e. A = 7ri? 4 ) of the heavy-light 
and heavy- heavy mesons for the non-supersymmetric background ( TTTT) . for m# = 1 
and mn — 7, with m^=0 in both cases. For large A, the heavy-heavy quark mass is 
suppressed. 



the 't Hooft coupling (depending on m#) above which the HL mesons are heavier than 
the HH mesons, which is unphysical from the point of view of QCD. Below the critical 
value the HH meson mass is larger than the HL meson mass. Accepting intermediate 
values of the 't Hooft coupling in the range of 1 < A < 50, it is possible to find an 
appropriate parameter region where realistic mass spectra are obtained. These spectra 
will be explored elsewhere. 

Moreover, we find numerically that for tq ~ R ~ A 1 / 4 large, the HL spectrum 
M 2 (R) approaches 2^ 7 ( u ' 2 — w i)\p=o, ie. is proportional to the distance of the two 
brane probes at p = in the far IR. Considering the embedding solutions in Fig. 4, we 
find that this mass is equivalent to the minimum energy of a string stretched between 
the two branes. At large 't Hooft coupling, this string is much shorter than a string 
stretched in between the same branes at p — > oo. Thus, the O(X ) contribution in the 
non-supersymmetric HL spectrum is much smaller than in the supersymmetric case, 
where it was proportional to mu — itll, see Eq. (T4"T|) . This implies in particular that 
in the non-supersymmetric case, the HL masses also tend to zero for A — > oo, though 
more slowly than the HH masses. 



5 Wilson loop for heavy-light mesons 

Unlike the spectrum of single-flavor mesons, the heavy-light spectrum of the super- 
symmetric theory does not vanish in the strong 't Hooft coupling limit, cf. Eq. (j4"8"j) . 
A better understanding of this behavior can be obtained by studying the forces (or the 
QCD-like tension) between the two quarks. For this, it is helpful to consider the quark 
anti-quark potential V q q which can be found by a standard Wilson loop computation 
similar to that in the case of single-flavor mesons [38J. 

The potential V q g is derived from the expectation value of a parallel Wilson- Polyakov 
loop, W = jjTrPe l -> A ° dt . In the dual gravity theory, it is represented as 

(W) ~ e~ s , (58) 
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E(L) 




Fig. 7: a) Numerical plots of the energy E(L) for HL mesons (A), LL mesons (B) and 
HH mesons (C). The circle at the endpoint of the curve (A) shows a finite string energy 
E at length L = 0. Here we set q = 5 and R — 1, and the brane positions are taken at 
r maxi = 10 and r max2 = 15, respectively, b) Schematic plot of the Wilson loop. 



with Nambu-Goto action 

S 



2na 



/drda^Z—dethab , (59) 



and induced metric h a b = G lxv d a X^d\ ) X 1 ' . The string world-sheet is parameterized by a, 
r, which in static gauge are set as X° —t — r and X 1 = x l = a. In the background (Q 
the Nambu-Goto Lagrangian becomes 



£ng = "&b / da e " /2A ^y r ' 2 + (^) 4 A2(r) ' (60) 

where the prime denotes the derivative with respect to a. We suppose here that the 
test string has a deformed U-shape whose endpoints are on the two D7 branes, as 
shown in Fig. [7b. 

Let us first consider the single-flavor case. From the Lagrangian (1601) . we find the 
relation 

e*/ 2 - 1 ( ^YA 3 (r) = h , (61) 

^{r/RYA 2 {r) + {dr/daf ^RJ 

where h denotes a constant of motion. Here we need to introduce two parameters, r min 
and r max : r min is determined by d a r\ rmin = and defines the bottom of the deformed 
U-shaped string, while r max is given by the position of the D7 brane, the endpoints 
of the string. Replacing h by r min by using the relation h = e^ 2 (-^) A 2 (r)\ rmin , we 
determine the energy E and the spatial distance L between the quark and anti-quark 
from 



Tmax j r^max 



L = 2R 2 I dr I L , E= — I dr I E , (62) 



A(r)e^' 2 



~ e^ mm ) r ± min A{r min y/ (e*WrM(r)4) ' 
1 

rM(r)yi^V 4 yl(r)V (e*^— ^ mm A{r mm Y) ' 
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(63) 
(64) 



where r min is restricted as < r min < r max . 

For the two-flavor case, L and E are given by 



Tmaxl f*T max2 



L = RU dr I L + drh), (65) 



mm 1 rmn 

fmaxl f > fmax2 



E = ^-;[ I <lrI L + I ,h 1 E ) . (GO) 



where r max \ and r maX 2 denote the positions of the light and heavy quark branes, re- 
spectively. They are equivalent to W\ and w-i for the supersymmetric case. 

For the supersymmetric background, numerical plots of the energy E(L) are shown 
in Fig. [Th. Curve (A) is a plot of E(L) for heavy-light mesons, while (B) and (C) 
correspond to those of the light-light and heavy-heavy mesons, respectively. The curves 
(B) and (C) end at (L,E) = (0,0); the energy is zero for vanishing quark- ant i- quark 
distance, E(0) = 0. However, as mentioned above, in the heavy-light case the energy 
remains finite, even for L = 0. 

For L = the string in Fig. [7b stretched between the D7 branes at r max i and r max2 
becomes just a straight line perpendicular to the boundary. Then, in the Nambu- 
Goto Lagrangian ([60]) . we write da = dr/r' and take r' = drjdo — > oo. In the 
supersymmetric case with q = we have A(r) = 1, $(r) = 0. Thus, from (1601 we get 



Tmax2 



Tmaxl 



This agrees with the meson mass result 

Note that the ordinary QCD string associated with the flux tube in between the 
quark anti-quark pair can be thought of as a projection of the Wilson loop on the 
boundary of the (asymptotic) AdS background, see Fig. [7b. The Wilson line for L = 
is aligned along the r-direction and therefore projected to a point on the boundary. 



6 Conclusion 

We have investigated heavy-light mesons in a holographic set-up by considering the 
non-Abelian Dirac-Born-Infeld action for two D7 brane probes embedded at different 
positions. The embedding matrix is chosen to be diagonal, whereas the heavy-light 
mesons arise from the off-diagonal elements of the fluctuation matrix. 

We considered both supersymmetric and non-supersymmetric backgrounds and 
found that the dominant contribution to the heavy-light meson mass is 0(1) in the 
't Hooft coupling A, whereas for the heavy-heavy and light-light mesons it is 0(l/\/~\). 
Although this result is unexpected from the point of view of QCD, it is consistent with 
our Wilson loop calculation for the heavy-light mesons, as well as with the effective 
field theory holographic approach to heavy-light mesons of [32J which uses the Polyakov 
string action. As discussed in section 4.2, the 0(1) contribution to the HL masses is 
proportional to the minimal energy of a string stretching in between two branes. In the 
non-supersymmetric case, the minimum energy corresponds to a string located close 
to the singularity (i.e. at p = in Fig. 4). At large 't Hooft coupling this energy is 
much smaller than m# — mi,- 
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As far as the spectrum of heavy-light mesons is concerned, we found that the HL 
spectrum lies below the HH spectrum for intermediate values of the 't Hooft coupling, 
1 < A < 50. At larger 't Hooft couplings, the HL meson masses exceed those of the 
HH mesons, which would be in conflict with phenomenology. Our results are consistent 
with scenarios in which QCD develops an infrared fixed point rather than a singularity, 
see e.q. [391 and references therein. In this greement with phenomenology would 

be achieved, if the gauge coupling a s {Q 2 ) were of order 1 or less at the fixed point 
(which corresponds to A = 4ira s (Q 2 )N c « 30). 

A new element of our approach as compared to [32] is that it allows to distinguish 
scalar from vector mesons. This opens up the possibility to study the heavy quark 
spin effect suppression known from QCD. We see the degeneracy of scalar and vector 
masses for large heavy quark mass in the M = 1 supersymmetric scenario. It will be 
interesting to study the non-supersymmetric case in the future. 

We conclude with some remarks on the range of validity of our approach of using 
the non-Abelian DBI action for two separated branes. Generally one may expect that 
the DBI is valid for branes separated at most by the string length l s . Here, however, 
our branes are separated by a larger distance if m# > 1 in our units of setting the string 
tension to one. Nevertheless our use of the non-Abelian DBI is justified by the fact 
that we obtain agreement with the classical string calculation of [32J, as well as with 
the semiclassical Wilson loop analysis of a string stretching between the two branes 
performed in section above in the present paper. Thus the classical analysis at larger 
length scales dominates over quantum fluctuations at shorter scales. This is at least 
in part due to the supersymmetry of the problem. We also note that restricting to 
heavy quark masses of order 0(oo) in the string tension will not affect the unusual A 
dependence, which intrinsically reflects the strong coupling behaviour. Physical quark 
masses correspond to m# <C 1. In this regime, the term m# -m^ causing the unusual 
A dependence is negligible. Moreover, for the non-supersymmetric approach it should 
be noted that the fact that the branes approach each other in the deep interior will 
not alter the A dependence either, since all values of the coordinate p contribute to the 
action. This is again in agreement with the results of [32] . 
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